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' Abstract 

We introduce in this paper a new interacting particle system intended to model an example of 
ecological succession involving two species: Pteridium aquilinum L. commonly called bracken 
and Fagus sylvattca L. or european beech. The species both compete over vacant sites, as- 
similated to windfall, according to fundamentally different evolutionary strategies. Our main 
objective is to exhibit phase transitions for this process by proving that there exist three 
\Q • possible evolutions of the system to distinct ecological balances. More precisely, the whole 

population may die out, the european beech may conquer the bracken, and coexistence of 
both species may occur. 
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1 Introduction 



Oh 

The ecological model that we introduce in this article is a Markov process in which the state at 
time t is given by some configuration rjt : Z d — > {0, 1, 2}. A site x S Z rf is said to be empty or 
vacant if T)(x) = 0, occupied by Pteridium aquilinum L. or bracken if r]{x) = 1, and occupied by 
Fagus sylvatica L. or european beech if 77(2;) = 2. If the site x is empty, it will be called sometimes 
a windfall. In what follows, we will formulate our results in the ecological context. However, to 
f^T) ' make the proofs more explicit, we will use another terminology, replacing the brackens by red 

particles and the european beeches by green particles. The evolutionary strategies of the two 
species involved in the competition are quite different and can be described as follows. 



1. The brackens are rhizomatical ferns that colonize their neighboring sites thanks to an under- 
ground plexus of rhizomes. Moreover, their optimal development requires sunlight so that they 
can colonize the windfalls only. 

2. The european beeches produce an important amount of seeds that are spatially dispersed by 
animals (zoochory) and can remain in the soil as seed bank ; the development occurs on the 
shadow only and so, close to the adult trees or sheltered from the brackens. 

With this in mind, we can now describe the dynamics of the ecological succession. Let V be the 
set of neighbors of in Z d given by V = { y E Z d , || y \\ < M } for some positive integer M and 
' some norm || • || on R d . 

1. A bracken (resp. a beech) standing at x tries to give birth to its progeniture on the neighborood 
x + V at rate Ai (resp. A2). The new individual is sent to a site chosen at random from the 
neighboring sites. 

2. If the target site is a windfall, the birth occurs. Otherwise, it is suppressed. 

3. A bracken gives way to a new beech at rate 1. 

4. Finally, a beech lives for an exponentially distributed amount of time with mean 1 and then 
dies out to give way to a windfall. 

For i € {1,2}, denote by Vi(x, rf) the proportion of neighbors of x occupied in the configuration 77 
by a type i particle. Then the points 1 and 2 are equivalent to the fact that a windfall becomes 
occupied by a particle of type i at rate A^ Vi(x, rf). 

We now formulate our results and construct step by step the phase space of the process in terms 
of the parameters Ai and A2. First of all, observe that if we start the evolution with a population 



American Mathematical Society 1991 subject classifications. 60K35 

Key words and phrases. Competition model, phase transition, interacting particle system. 



2 



N. LANCHIER 



of european beeches in the absence of brackens, the process we obtain is called the basic contact 
process with parameter A 2 . We know that, in this case, there exists a critical value A c G (0, oo) 
such as the following holds. If A 2 < A c then the process converges in distribution to the all empty 
state, else there exists a stationary measure that concentrates on configurations with infinitely 
many beeches. See for instance Liggett (1999), Sect. 1.2. 

Intuitively, we can think that, on the set where coexistence of both species does not occur, this 
result goes on holding for the ecological succession model. If A 2 < A c , it is obvious, i.e. both 
species die out or coexist. If we now suppose that Ai > A 2 > A c , we can prove that starting 
with infinitely many occupied sites, the process converges in distribution to a stationary measure 
that concentrates on configurations with infinitely many beeches. To do this, consider two initial 
configurations r/ 1 and r\ 2 such as r\ 2 (x) = if r/ 1 (x) = 0, and r\ 2 (x) = 2 otherwise, i.e. r\ 2 can 
be deduced from n 1 by replacing the brackens by beeches. Observe that this implies that % is 
the basic contact process with parameter A 2 . Then by using a standard argument due to Harris 
(1972) we may run both processes simultaneously in such a way that if % (x) = 2 then r\l (x) ^ 0. 
Hence, we can finally state the following 

Theorem 1 Starting with infinitely many brackens, the process r\ t converges in distribution 

a. either to the all empty state or to a measure that concentrates on configurations with infinitely 
many brackens (and beeches) if A 2 < A C7 

b. either to the upper invariant measure of the supercritical contact process with parameter 
A 2 or to a measure that concentrates on configurations with infinitely many brackens (and 
beeches) if M > A 2 > A c . 

Nevertheless, such a coupling fails to prove an analogous result in the case Ai < A 2 . But, since 
the particular dynamics of the ecological succession is inclined to favour the development of the 
beeches, it looks reasonable to think that, in the case A 2 > A c , the beeches take over provided that 
Ai < A 2 . In fact, relying on some technicalities introduced by Neuhauser (1992), we can prove the 
two following theorems. 

Theorem 2 We suppose that A 2 > A c and i] is translation invariant. 

a. If Ai < A 2 then, starting the evolution with infinitely many beeches, the brackens die out, 
i.e. i] t => H2 the upper invariant measure of the contact process with parameter A 2 . 

b. If Ai = A 2 then, starting the evolution with infinitely many occupied sites, r\t => M2- In 
particular, the beeches take over. 

The last step is to show that the ecological succession does not reduce to the basic contact process, 
i.e. that coexistence of both species occurs for an open set of the parameters Ai and A 2 . If A 2 = 
and d = 2, the process is known as the forest fire model. See for instance Durrett and Neuhauser 
(1991). In this context, = alive, 1 = on fire and 2 = burnt. We start the evolution with one 
burning tree in a virgin pine forest and follow the spread of the fire. A burning tree sends out 
sparks at rate Ai = a to one of its nearest neighbors x + ei, x + e 2 , x — e\ or x — e 2 and burns 
for an exponentially distributed amount of time with parameter 1. Finally, a new pine regrows 
spontaneously on the burnt site at rate (5 > 0. By comparing this particle system with a one- 
dependent oriented percolation process, Durrett and Neuhauser (1991) have proved that, for fixed 
f3 > 0, there exists a critical value a c € (0, oo) such as the fire has positive probability of not going 
out if a > a c , i.e. there exists a nontrivial stationary distribution. Relying on their result and 
using an idea of Schinazi (2001) wc can prove that this goes on holding for the ecological model 
with small values of A 2 > 0. More precisely, we have the following 

Theorem 3 We suppose that d = 2 and that V reduces to the nearest neighbors. If Ai > a c there 
exists a critical value (3 C € (0, oo) such as coexistence occurs as soon as A 2 < (3 C , i.e. there exists 
a stationary distribution that concentrates on configurations with infinitely many brackens. 
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In conclusion, there exist three possible evolutions of the system to distinct ecological balances. 
That is, depending on the parameters Ai and A2, the whole population may die out, the beeches 
may win the competition, and coexistence of both species may occur. For an picture of the space 
phase, see Fig. 1. On the other hand, it looks reasonable to conjecture, in view of the ecological 
dynamics, that for any Ai > a c 

1. There exists a critical value (3 C (Ai) G (0, 00) such as conclusion of Theorem 3 holds if and 
only if A 2 </3 c (Ai). 

2. The critical value (3 C (Ai) is strictly increasing with respect to Ai. 




empty 



Fig. 1 



Putting things together, the phase space of the particle system can be divided into three conncx 
components corresponding to the three possible evolutions of the process. See Fig. 2 for a pic- 
ture. This also implies that there exists a critical interval I c such as for any Ai G I c , the three 
possible evolutions may occur. In particular, we can notice that for any Ai G I Cl an increase of the 
beeches birth rate A2 can make disappear the brackens and the beeches themselves. In the case 
d = 1 and V = { — 1, 1}, we have observed this phenomenon by simulation for any Ai > a c with 
a c ~ 7, 95 x card V = 15, 9. In particular, in view of the special geometry of the one dimensional 
process, we think that coexistence cannot occur in d — 1 whenever A2 is supercritical. 

The article is organized as follows: In Sect. 2, we describe in greater details the graphical con- 
struction of the dual process using an idea of Harris (1972). In Sect. 3, we prove Theorem 2 for 
Ai < A 2 in any dimension relying on the graphical construction given above and some results of 
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Neuhauser (1992). The proof of Theorem 2 for Ai = A2, essentially based on classical properties 
of random walks on Z d , is carried out in two steps. In view of the well-known specificities of 
symetric random walks depending on the dimension, we start, in Sect. 4, by proving the theorem 
in d < 2, and conclude, in Sect. 5, dealing with the case d > 3. Finally, Sect. 6 is devoted to the 
description of results of Durrett and Neuhauser (1991) and the proof of Theorem 3. 




Fig. 2 



2 Construction and properties of the dual process 

We begin by constructing the process from a collection of Poisson processes in the case Ai < A 2 . 
For x, y £ Z d , x — y E V, let { T* ,v ; n > 1 } and { U% ; n > 1 } be the arrival times of Poisson 
processes with rates A2 • (cardV) -1 and 1 respectively. At times T*' y , we draw an arrow from x 
to y, toss a coin with success probability (A 2 — Ai)/A 2 , and, if there is a success, label the arrow 
with a 2. If a; is occupied, y is vacant and the arrow is unlabelled, then a birth will occur in y. If 
the arrow is labelled with a 2, then the beeches only can give birth through this arrow. Finally, 
at times U%, we put a x at x. In view of the dynamics, the effect of a x is to kill beeches and 
to permute the brackens into beeches. A result of Harris (1972) implies that such a graphical 
representation can be used to construct the process starting from any 770 <G {0, 1, 2} z . 

We say that [x, 0) and (y, t) are strongly connected, which we note [x, 0) — > {y,t), if there is 
a sequence of times So = < s\ < s 2 < ■ ■ ■ < s n < s„+i = t and spatial locations xq = 
x, xi, X2, ■ ■ ■ x n = y so that 
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(1) for 1 < i < n, there is an arrow from Xi-\ to Xi at time and 

(2) for < i < n, the vertical segments {x^ x (si, s i+ i) do not contain any x . 
If instead of (2) we have the condition 

n 

(3) the set {xi} x (s i; s i+ i) contains exactly one x, 

i=0 

we say that (x, 0) and (y,i) are weakly connected, which we note (x, 0) — (y, t). 

After constructing the graphical representation we now define the dual process. Since the x do 
not kill the brackens but convert them into beeches, we must take into account the paths that 
contain at most one x. So, to construct the dual process, we reverse the arrows and time by 
mapping s = t — s, and let 

f] x t = {yeZ d ; (x, 0) - (y, t) or (x, 0) - (y, t) }. 

Since it is in general easier to work with a forward process, we replace fjf by the dual process fff 
that is constructed by restoring the direction of time and arrows and by letting 

ff t = { y e Z d ; (x, 0) -> (y, t) or (x, 0) - (y, t) }. 

As in the basic contact process, fjf and f/f have the same distribution. 

Before sketching the proof of Theorem 2, which is based on the graphical representation, we give 
some definitions intended to describe the geometry of the dual process. First of all, note that 
{f)f;0<s<i} has a tree structure composed of the points that are either strongly or weakly 
connected with (x, 0). To figure on this distinction, denote by L the tree formed by the set of 
points in Z d x R + that are strongly connected with (x, 0), i.e. 

L = {(y,t)eZ d xR+; (x, 0) -> (y, t) }. 

Clearly, the growth of L is broken at some points by a x from which a new tree takes form so 
that we obtain the picture of a pyramid of trees connected with L by some x . In the following, 
r will be called the upper tree starting at (x, 0) and the trees starting at a x the lower trees. For 
an illustration, sec Fig. 3. To motivate such a breaking up, note that, in view of the translation 
invariance of the graphical representation in space and time, the upper tree and the lower trees are 
identically distributed. Furthermore, T has obviously the same distribution as the tree structure of 
the contact process with parameter A2, that is the set of points in Z d x R + that can be connected 
by some path to (x,0). See for instance Durrett (1995), Sect. 3 or Liggett (1999), Sect. 1.1 for 
a construction of the basic contact process. So, one of the main interest of our decomposition is 
that it allows to break up the tree starting at (x, 0) into identically distributed pieces that have 
some usefull similarities with some well-known structures. 

Now, denote by rfc and fit' 2 the subsets of the dual process given by 

77^ = {yeZ d ; (z,0)->(y,t)} and ft 2 = { y e Z d ; (x, 0) - (y, t) }. 

By analogy with the breaking up of the tree structure in upper tree and lower trees, the elements 
of fjf' 1 and fit' 2 will be respectively called upper ancestors and lower ancestors. As in the mul- 
titype contact process, the tree structure of the dual process {r/f ; < s < t} allows to define 
an ancestor hierarchy in which the members are arranged according to the order they determine 
the color of (x,0). Here, the special geometry of the dual described above plays a leading role 
since the color of (x, 0) also depends on the nature of the ancestors. To specify this idea, we now 
describe in greater detail the hierarchy. 
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First of all, suppose that Ai = A 2 . Denote by fj^(n) the n-th member of the ordered ancestor set 



and let (k) = rj? (nk) be the k-th upper ancestor. Later on, (1) will be called the distin- 
guished particle, and the first ancestor the first particle. If the first particle is a upper ancestor 
that lands on a red particle (resp. a green particle), then it will paint {x,0) in red (resp. green). 
If it is a lower ancestor that lands on a red site, in view of the passage through one x , it will bring 
a green particle to (x, 0). If it is a lower ancestor that lands on a green site, the path taken by the 
green particle until it is destroyed by a x can prevent the rising of some ancestors toward (x, 0). 
Hence, we discard all the ancestors whose path crosses the green particle. Then, we repeat the 
same procedure with next remaining ancestor, and so on. 



M) 



M) 



1 2 [_n 1 

2 16 5 



Fig. 3 



Fig. 4 



If A2 > Ai, the arrows labelled with a 2 are forbidden for the red particles. So, if fjf (1) docs not 
cross any 2-arrows or lands on a green site or a windfall, we apply the same procedure as before. 
If it is a lower ancestor that lands on a red site and that crosses a x before the first 2-arrow, it 
will paint (x, 0) in green. In the other cases, we follow the path taken by the particle until we first 
cross a 2-arrow. Then we discard all the ancestors of the point where this arrow is attached, and 
start afresh with the first remaining ancestor, and so on. For instance, in Fig. 4, the first particle 
and the third ancestor are lower ancestors that land on a green site so they cannot paint (x, 0) 
any color. Moreover, the path of the second ancestor is forbidden by the rising of a green particle 
at site x + 1 so that the distinguished particle i)f (4) paints (x, 0) in red. If we now suppose that 
A 2 > Ai, wc can notice that r?f (4) also fails and that the fifth ancestor paints (x, 0) in green. 

3 Proof of Theorem 2 for A : < A 2 

In this section, we will prove Theorem 2 for Ai < A 2 in any dimension. The strategy of the proof is 
quite simple ; it relies on a usefull idea of Neuhauser (1992) that consists in breaking up the path 
of the distinguished particle at certain points into i.i.d. pieces. We will first remind this result, 
describing the intuitive idea on which it is based, and then conclude that the brackens die out as 
soon as Ai < A 2 . 



To begin with, observe that if the upper tree starting at (x, 0) does not live forever then, for t 
large enough, rjf — i)*' . In particular, in view of the passage of any lower ancestor through one 
x, (x, 0) cannot be reached by a red particle (see the description of the ancestor hierarchy in Sect. 
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2). So we can suppose from now on that the upper tree T lives forever. Note that the probability 
of such an event is given by the survival probability of the contact process with parameter A 2 
starting from one infected site, that is obviously positive since A 2 is supercritical. In this case, we 
proceed in two steps by firstly checking that the path of the distinguished particle is forbidden 
for the red, and then proving that we can bring a green particle to (x, 0). To do this, we define 
an embedded random walk for the distinguished particle by breaking its evolution at some points 
called renewal points. 

We now specify this construction. First of all, denote by Sk, k > 1, the jumping times of the 
distinguished particle. Then, follow the path of the particle inside the upper tree and, whenever it 
jumps to a new site ^^(l), look at the branch starting at (^^(l), s/c). If this branch lives forever, 
WsfcHl)' s fc) wm be called a renewal point. Let (S° , T°) be the location of the n-th renewal. For 
an illustration, see Fig. 5. Denote by X° the spatial displacement between consecutive renewal 
points, and by t° the corresponding temporal displacement so that 

n n 

S n = x + X * aild T « = T *°- 

j=l i=l 

Our main ingredient to prove Theorem 2 is given by the following proposition. We just describe 
the intuition behind this result. For the details of the proof, see Neuhauser (1992), Sect. 2. 

Proposition 3.1 If the upper tree lives forever, {(X° , t° )}i>i form an i.i.d. family of random 
vectors on Z d x R+. Moreover, the tail distributions of X° and t° have exponential bounds, i.e. 
there exist positive constants C and (3 such as 

P(\X°\>t) < Ce- pt and P (t° > t) < C e~ 0t 

holds for all t>0. 

Sketch of the PROOF. First of all, denote by Co the first jumping time of the distinguished 
particle, that is cto = inf { t > ; rfc' (1) hits a x }, by x\ its spatial location after <ro, and by (3\ 
the branch of T starting at (xi,o- n ). For a picture, see Fig. 6. If (3\ lives forever then (xi,cr n ) is 
the first renewal point. Else, we define the sequences {(xk, Ck-i)}k>i and ((3k)k>i as follows. If 
(7fc_i < oo, let Xk be the location of the particle after o~k-i and (3k the branch starting at (xk, <Tfc-i). 
Note that such a branch always exists since we have supposed that T lives forever. Then denote 
by o~k the time when (3k dies out. The sequences are defined until &k is equal to infinity. If o~k = oo 
then (3k lives forever and {xk, e>fc_i) is the first renewal point of the distinguished particle. To find 
the next one, we start over again the whole procedure replacing (x,0) by (xk,ak-i), and so on. 
Now, look at the path of the particle between time and time o~k-i. Whenever the particle jumps 
within a new branch Pi, the location of its target site Xi clearly depends on the past history of T, 
that is on the geometry of T between and <7i_i. On the contrary, from time <7fc_i, the particle 
stays forever inside (3k and its path only depends on this branch so that what happens before and 
after a renewal point is determined by disjoint parts of T. This implies that the random vectors 
(X° , Tj° ) are independent. Finally, since the graphical gadget is translation invariant in space and 
time the vectors (X° ,t°) are also identically distributed. The proof of the exponential bounds, 
more technical, can be found in Neuhauser (1992), Sect. 2. □ 

By Proposition 3.1, the sequence (S°,T°) defines an embedded random walk for the particle. 
Moreover, the exponential bounds on the displacements X° and t° give us control over the loca- 
tion of the particle between consecutive renewal points. Since the branch starting at each renewal 
grows at most linearly in time, we know that the particle stays within a triangle whose base con- 
tains the following renewal. For a picture, see Fig. 8. 



8 



N. LANCHIER 



A triangle is said closed for the red if the distinguished particle goes through an arrow labelled 
with a 2 to cross it. Since Ai < A 2 , the event that a triangle is closed for the red obviously occurs 
with positive probability. So, by the Borel-Cantelli Lemma, choosing t large enough, we can find 
arbitrarily many closed triangles. 



(x,0) 



2 



(SIX) 9 



(SlT°) o 



<?1 



02 



<7:\ 
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{SIX) ? * 



Xi x 2 x 3 



X4 



Fig. 5 



Fig. 6 



To conclude the proof, we now use the dual process fj/'* , < s < t, starting at (x,t) and 
determine the ancestor hierarchy after t units of time by going backwards in time. We construct 
by induction a sequence of upper ancestors (k), k > 1, that arc candidates to bring a green 
particle to (x,t). The first member of the sequence is the distinguished particle ^^(l). Then, 
we follow the path the distinguished particle takes to paint (x, t) by going forward in time until 
we first cross a 2-arrow. Finally, we look backwards in time starting from the location where this 
arrow is attached and discard all the offspring of this point (that is the next few members of the 
ordered ancestor set) and all the ancestors that land at time on the same site as the distinguished 
particle. The second member of the sequence is then the first upper ancestor that is left. We repeat 
the procedure with this ancestor, and so on until we run out of upper ancestors. Now, denote by 
£t the set of members of the sequence and by B s "' t} the set of sites occupied at time s by a green 
particle. Since the tree T is linearly growing in time, we can make the cardinality of £t arbitrarily 
large by choosing t large enough. In particular, given e > and M > there exists a time to > 
so that P ( card (£t) < M ) < e for any t > t^. By Lemma 9.14 of Harris (1976), this together with 
the translation invariance of tjq and the fact that r)o(x) = 2 with positive probability, imply that 

lim P(6-i n £? 1 <x ' t) = 0) = 0. 

Hence, for t large enough, at least one candidate lands on a green site. Now, consider the first 
member of the sequence ^^(k) that lands on a site occupied by a green particle and follow the 
path of this particle by going forward in time. Then, whenever it crosses an arrow, its target site 
is either empty or already occupied by a green particle (resulting from the passage of a red particle 
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through the first x located under the arrowhead). Thus, (a;, t) will be eventually painted in green. 
This completes the proof of Theorem 2 for Ai < A 2 . 

4 Proof of Theorem 2 for Ai = A2 in dimension < 2 

The technicalities and tools we will make use to prove Theorem 2 in the case Ai = A2 are quite 
different depending on the dimension of the state space. In this section, we will deal with the case 
d < 2, relying on the recurrence of 1 and 2 dimensional random walks. To begin with, we will show 
that the first particle can be trapped with probability one inside a lower tree that lives forever, so 
that, for t large enough, the first ancestor is a lower ancestor. At this point, the worst scenario 
we have in mind is that the distinguished particle lands on a red site and the first ancestor on a 
green one. In such a case the first particle, which is a lower ancestor, cannot paint (x, 0) any color 
whereas the distinguished particle can possibly bring a red particle to (x,0). To conclude, we 
will then prove that this bad event is negligible showing that, with probability one, we can make 
coalesce the distinguished particle and the first particle and so make them land on the same site. 
If they both land on a green (resp. on a red), the distinguished particle (resp. the first particle) 
will paint (x, 0) in green, and the theorem will follow. 

From now on, we suppose that Ai = A 2 and denote by A their common value. The first step is 
given by the lemma 

Lemma 4.1 There exists a time 61 a.s. finite such as for any t > 6\, the first particle is a lower 
ancestor. 

Proof. To begin with, let Sfe, k > 1, be the jumping times of the first particle and Xk its location 
before Sk- Denote by o\ the first time the particle crosses a x, i.e. 

o\ = inf { t > ; fjf (1) is a lower ancestor }, 

and by fii the lower tree starting at (xi,t7i), that is the first lower tree the particle visits. See 
Fig. 7 for an illustration. In view of the tree structure, once the particle penetrates in fii, it 
remains trapped inside (as Ions as the lower tree is alive). Hence, if fii lives forever, 17^ (1) is a 
lower ancestor for any t > a\, and the proof is done. Else, denote by a 2 the first time the first 
particle visits a new lower tree after Sli dies and by Q.2 the lower tree starting at (17^ (1), 02). 
Note that for all k > 1 the path the particle takes to climb from (xk,s^) to (x,0) contains one 
x so (72 is finite a.s. and fi2 is well defined. While the particle is not trapped in a lower tree 
that lives forever, we thus construct by induction a sequence of trees fi& visited by rjf (1). Now, 
denote by B n the event that the first n trees fii, Q2, ■ ■ ■ , fin are bounded and, for any k > 1, by 
Ak the event that the fc-th tree lives forever. If Ak does not occur then fifc+i is well defined and 
the event Ak+i is determined by parts of the graph that are after fi/j dies so Ak and Ak+i are 
independent. More generally, since the trees fii, fi 2 , • • • , fifc+i are disjoint, A\, A2, . . . , Ak+i are 
independent. Moreover, fi^ has the same distribution as the tree structure of the contact process 
so the probability that Ak occurs is given by p\, the survival probability of the contact process 
with parameter A starting from one infected site. This implies that 

n 

P(B n ) =F(i;n...n^n^) = *[[p{AI) = (i-paT- 

fe=i 

Finally, since A is supercritical, the survival probability of the contact process p\ is strictly positive 
so that lim P (B n ) — 0. This completes the proof of the lemma. □ 

n— >oo 

The next step of the proof is to show that the particles coalesce with probability one. To do this, 
note that after penetrating in a lower tree Ti that never dies, the first particle is weakly connected 
with (x,0) and so jumps within a new branch each time it meets a x. In particular, from time 
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01, the path of the particle can be broken up, as for the distinguished particle, into i.i.d. pieces. 
We define the renewal points of the first particle as before replacing the upper tree T by the lower 
tree Ti starting at (Sq , 9\). Here, Sq is the site where the particle jumps at time 6\. We denote 
by (S^ , ) the location of the n-th renewal after 6\ , and by X- and r/ the spatial and temporal 
displacements between two consecutive renewals, so 



E 

i=l 



X 



and 



T = 



i=l 



By translation invariance of the graphical representation (see the description of the tree struc- 
ture in Sect. 2), the families {(X° , t° )}i>i and {(X^ , )}i>i are identically distributed so that 
Proposition 3.1 remains valid for the random vectors (X^ ,t^). In particular, as long as their 
triangles do not collide, both particles behave nearly like independent random walks. This con- 
stitutes the main ingredient to prove coalescence. We now specify this idea in greater details, 
describing the approach of Neuhauser (1992), Sect. 3. 
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The first idea is to extend the notion of renewals for both particles, that is to break up the set of 
both paths into i.i.d. pieces. To do this, we say that an ancestor is good at time t if it did not 
meet any arrow since its last renewal. Observe that if both particles are good at the same time, 
what happens before and after that time uses disjoint parts of the graph and then is independent. 
We now prove that both particles are good i.o. at the same time. 

Lemma 4.2 P ( both particles are good at the same time i.o. ) = 1. 

Proof. To begin with, we construct by induction two sequences of subscripts (nk)k>o an d {rrik)k>i 
as follows. We let no = 1, and for any k > 1 



m k = min { m > 1 , T r ° > T„ 1 fc _ i } 



and 



n k = min{n> 1 , T„ x > T° }. 



Then, denote by A k the event that the first particle lives without giving birth between time 2 n ra 1 fc _ ] 
and T^ k . Note that if A k occurs then both particles are obviously good at time T^ k . Moreover, 
in view of the exponential bound given by Proposition 3.1, for any T > 

PiT^-T^yT) < P ( T^ k - T^ lk _ 1 > T ) = P(r 1 °>T) < C e^ T 

so that P{A k ) > ( 1 - C e~ liT ) e - (1+ ^ T . Here, e~( 1+A ) T is the probability that the first particle 
lives without giving birth for T units of time. Since this holds for all T > 0, there exists a constant 
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e > such as P(Ak) > £■ On the other hand, the events Ak are determined by disjoint parts 
of the graphical gadget so they are independent. Finally, by the Borel-Cantelli Lemma, we can 
conclude that 

P ( both particles are good at the same time i.o. ) > P ( lim sup A/j ) = 1. 
This proves the lemma. □ 

To make coalesce both particles together, we now proceed in two steps. First of all, relying on the 
recurrence of 1 and 2 dimensional random walks, Neuhauser has proved that with positive prob- 
ability we can bring both particles within a finite distance K without collision of their triangles. 
See Neuhauser (1992), Sect. 3, for a proof. Then, as soon as the particles are close enough to each 
other, we try to make them coalesce. More precisely, we have the following 

Lemma 4.3 Let K be a positive integer and suppose that both particles are within a distance K 
at some time t>6\. Then the event B that the particles coalesce has positive probability. 

Proof. To find a lower bound for P(B), we first require the distinguished particle to keep still for 
3dK units of time. Then we can make both particles coalesce in less than dK steps by increasing 
or decreasing each coordinate of the first particle. To estimate this event, note that the probability 
of having neither birth nor death between time t and t + 1, a good oriented arrow between time 
t + 1 and t + 2, and a death between time t + 2 and t + 3 is given by 

e- 2 (1-e- 1 ) e- x {V^ 1 (l-e- x ). 

Since it takes the first particle at most 3dK units of time to reach the distinguished particle, we 
can conclude that 

P(B) > e -3^(i+A) [ e -2 (i- e -i) e ~ x \V\^ (l- e - x )] dK > 

where e - ?,dK ( 1 + x ) j s the probability that the distinguished particle survives without giving birth 
for 3dK units of time. This completes the proof. □ 

If we do not succeed in gluing the particles together, we use the restart argument given by Lemma 
4.2, i.e. we wait until both particles are good at the same time and then start over again the 
whole procedure. Since the set of paths is broken into i.i.d. pieces, we can apply the Borel- 
Cantelli Lemma to conclude that coalescence eventually occurs. 

We prove Theorem 2 in d < 2 by using, as before, the dual process fj/'* , < s < t, starting at 
(x,t). Since both particles coalesce a.s., we can suppose, taking t large enough, that they land at 
time on the same site. If it is a red site, the first particle, that is a lower ancestor by Lemma 4.1, 
will paint (x,t) in green. If both particles land on a green site, we have already proved in Sect. 
3 that the distinguished particle can bring a green particle up to (x,t). Finally, if the target site 
is a windfall, we start over again with the second ancestor, and so on. Since the tree starting at 
the coalescence point is linearly growing in time, we eventually find, by Harris (1976), an ancestor 
that lands on an occupied site and paints (x, t) in green. 

5 Proof of Theorem 2 for Ai = A 2 in dimension > 3 

The strategy of the proof to deal with the case d > 3 is quite different. To begin with, we will 
construct by induction an ordered set of ancestors, 7f (fc), k > 1, that are candidates to paint 
(x,0) in green. Using Lemma 4.1, we will prove that for any k > 1, and for t large enough, (k) 
is a lower ancestor that comes before the distinguished particle in the ancestor hierarchy. Then, 
relying on the transience of d-dimcnsional random walks for d > 3, we will extract a subsequence 
of particles jf (fcj), i > 1, that never coalesce. In particular, the number of sites occupied by the 
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particles can be made arbitrarily large so that we can eventually find a particle landing on a red 
site that will paint (x, 0) in green. 

We start by constructing inductively the particle system (fc), k > 1. The first member of the 
sequence is the first particle. Then, follow its path until it penetrates in a lower tree Fi that lives 
forever. By Lemma 4.1, we know that this occurs at some time 9\ a.s. finite. Now, consider 
the ancestor hierarchy at that time and discard all the ancestors that land on S a or on a branch 
that does not live forever. We remind that S is the spatial location of the first particle at time 
01. The second member of the sequence is then the first ancestor that is left. Observe that such 
an ancestor obviously exists since we have supposed that T lives forever. Moreover, using the 
arguments of 4.1, we can prove that it penetrates a.s. in a lower tree 1?2 that never dies. We then 
repeat the same procedure to define the third particle, and so on. For any k > 1, let F& be the 
infinite lower tree visited by the fc-th particle, 9k the first time the particle penetrates in and 
Sq the spatial location of the particle at time 9k- As for the first particle, we can break up the 
path of the k-th particle inside T k into i.i.d. pieces. We denote by (S£,T*) the n-th renewal 
point, by X* the spatial displacement between consecutive renewals, and by r 4 * the corresponding 
temporal displacement, i.e. 

n n 

^ = So + E X * T n = 9 * + E T ^ 

i=l i=l 

The next step is to show that for t large enough we can find arbitrarily many sites occupied by 
the particles (k), k > 1, i.e. there are infinitely many particles that do never coalesce. Before 
giving the idea of the proof, we start by proving some preliminary results intending to describe 
the behaviour of the particles. We give the proof for the first particle only but the same holds for 
the other members of the system. 

Lemma 5.1 Let K be a positive integer and Bk the spatial box [— K, K] d . Then for all t > 9i 
and z G Z d the event A = { the first particle leaves z + Bk in less than one unit of time and then 
survives without giving birth until t + K } has positive probability. 

Proof. If the particle stands out of z + Bk at time t, then P(A) can be bounded from below by 
e -(A+i)Jf {.^at is the probability that the particle lives without giving birth for K units of time. 
Otherwise, we require the first particle to increase or decrease its first spatial coordinate until 
leaving z + Bk- Taking 8 > such as (K + 2) 6 < 1, a straighforward calculation shows that 

P(A) > e- {1+x ^ K [e- 25 {l-e- 5 ) e- xs \V\- X {l-e- x5 )] K > 

since it takes at most K steps for the particle to reach a face. This concludes the proof. □ 



Lemma 5.2 For t > 6i and K > 0, denote by H ty K the event that the triangles of the first particle 
are contained in some box z + Bk from time t and for K units of time. Then, for any e > 0, there 
exists a large enough K so that Ht,K occurs with probability at least 1 — e. 

Proof. To make the notations easier, we will omit, all along the proof, the superscripts 1 that 
refer to the first particle. Moreover, since the graphical representation is translation invariant, we 
can suppose that the first renewal point after time t stands on the origin and let z = 0. Observe 
that in this case, H t ,K occurs if the spatial locations of the particle at the renewals do not leave 
the box i Bk and if each triangle is smaller than i K. Denoting by -% the location of the random 
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walk S n at time s and by N tl K the number of triangles between time t and t + K, we then obtain 

P(H tK does not occur ) < P ( N t K > — K 

\ ' m 

( K 2 

+ P the largest triangle is > — ; N t k < — K 
\ 2 m 

+ P ^ l7.IL > y for some t < s < t + K ; N t , K <^K 

where m = En. To begin with, we can bound the first term on the right-hand side by using the 
large deviation estimate 



P (^N tiK >^-K^j < Ce~ pK . 



Moreover, since the random vectors Xi, i > 1, have the same distribution, the second term can 
be bounded as follows 

/ K 2 

P the largest triangle is > — ; N t k < — K 
\ 2 m 

< p( max \\XA > -) < - K P ( \\XA > -) < -KCe~^ K . 

Observe now that if the random walk S n leaves the box | Bk then at least one of its coordinates 
is bigger than \ K so the third term can be bounded by 

rf / K 2 

P ( l^ <l> I > ~2 for SOmC 1 ~ s ~ 1 + K ' Nt ' K - ~ K 

i=l ^ 

(K 2 
|7 S (1) | > — for some t < s < t + K ; N t K < — K 
2 m 

where the superscript i refers to the i-th coordinate. We can use the reflection principle to bound 
this last term by 

On the other hand, Chebyshev's inequality gives for any 9 > 



(rs \ *™ n m 



2m,- 1 K 
i=l 



where <p(8) is the moment generating function of X^ 1 ' . Since EIj 1 ' 1 = and Var X^ < oo we 
can state that log (f)(6) < C 9 2 for some C > and for 9 small enough. In particular, taking 

9 = — ^= in the last expression we conclude that 

p ( Wlsi > — for somc t < s < t + K ; N tK < — K) < 2de~i VF 
\ °° 2 m J 

for K large enough. Putting things together, we can finally maintain that 

P (H t k ) > 1 - Ce~ pK - —KCe-^ K - 2de~^. 

m 

This completes the proof of the lemma. □ 
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We now prove that for t large enough the number of sites occupied at time t by the particles 7f (fc) 
can be made arbitrarily large. The strategy of the proof consists in extracting a subsequence of 
ancestors that do never coalesce. To do this, we proceed by induction. Let n > 1 be a positive 
integer and suppose that the n particles (fc,), 1 < i < n, never collide. Our main ingredient to 
prove that there exists a.s. another particle 7f (£) that does not coalesce with the first n ones is 
given by the following lemma. 

Lemma 5.3 For any e > 0, there exists a large enough K n such as the following holds: If for all 
1 < i < n both particles 7f (fci) and jf(£) are good and get separated by a distance > K n at some 
time t>0 then 

P ( the particles do never coalesce ) > 1 — s. 

Moreover, on the set where coalescence does not occur, the distance between the particles tends to 
infinity as t — > oo. 

Proof. Let 1 < i < n and suppose that the particles 7f (h) and -ff (£) are good and separated 
by a distance > K at some time t > 0. Then, Lemma 5.5 of Neuhauser (1992) implies that there 
exists a constant C > such as 

P ^the particles *y£ (fcj) and ~f£(£) do never coalesce) 

= P( lim Ik? (ki) - 7? = +00) > 1 - CK~ 1/la - 2CK~ 3/3 \ 

In particular, for K n large enough, 

P (^it(ki) and Jt(£) coalesce for some 1 < % < n^j 

n 

< p {lt{h) and -r?(t) coalesce) < nC K n 1/U ' + 2nC K n 3/32 < e. 

i=l 

This proves the lemma. □ 

Putting the last three lemmas together, we can prove the following 

Lemma 5.4 Let A e ^ K = { for every 1 < i < n, both particles jfih) an d 7?W ar e good and get 
separated by a distance > K at least once between time 9g and }. Then for K large enough, 

the event Ai^k occurs for infinitely many £>1. 

Proof. Since the first n particles do not coalesce and that lim^ = +00, we can find by Lemma 
5.3 a large enough £ such as the distances between the particles at time 9i are bigger than 4^/dK. 
Then, denote by Hi the event that the triangles of the ki-th particle are contained, between time 
9i and 6e + K, in some box fli = Zi + Bk- Observe that £ has been chosen so that VLi fl Clj = 
as soon as i 7^ j. Hence the events Hi are determined by disjoint parts of the graph, and then are 
independent. This together with Lemma 5.2 implies that for K large enough 

P {H l n h 2 n • • • n H n ) > (i-e) n . 

Now that the first n particles are trapped inside large disjoint cubes, we require each of them to 
be good at least once between time 9i and 9t + K. Since this occurs if the n particles have at least 
one renewal in this time laps, the probability of this event can be bounded from below by 

n 

\\ P{t" 1 < K- 1) > [l-Cexp(-/3(Jr-l)) 

i=l 

where C and (3 are some positive constants. The last thing we need is that the £-th particle is 
good and its distance from each other ancestor is > K between 9i + 1 and 9t + K. Since the 
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cubes 0, are separated by a distance > 2K , this obviously occurs, by Lemma 5.1, with positive 
probability. Putting things together we can state that there exists an r\ > such as P(Ae^) > V 
for I large enough. Observing finally that the events Ag lt K and A^ 2i k are independent as soon 
as \0£ 1 — 0£ 2 \ > K, we can conclude by the Borel-Cantelli Lemma that a.s. the events Ai^k occur 
for infinitely many I > 1. □ 

On the event Ag^x coalescence occurs with probability at most e for K large enough. Moreover, 
Ae t K occurs for infinitely many I > 1, so we eventually find a particle ^(i) that never coalesces 
with the first n members of the sequence. Since n has been chosen arbitrarily, this proves that 
there arc infinitely many particles that do never coalesce. To conclude the proof of Theorem 2, we 
use again the dual process fj/'* , < s < t, starting at (x, t). Since, by choosing t large enough, we 
can find arbitrarily many different sites occupied by the particles "f£(k) and that rj Q is translation 
invariant, there exists a.s. a lower ancestor 7f (n) that lands, at time 0, on a red site. Now, look 
at the ancestors that come before -ff (2) in the hierarchy until we find a particle that lands on an 
occupied site. If it is a red site, it will bring a green particle to (x, t) and the proof is done. If it is 
a green site, Sq will be vacant at time (t — 9\) + whatever the color of the next ancestors. Then, 
we look at (2), and so on. If however not one of the ancestors that come before 7f (n) succeeds 
in painting (x, t) in green, this last one will do it. This completes the proof of Theorem 2. 

6 Proof of Theorem 3 

This last section is devoted to the proof of Theorem 3. In particular, we will prove that, in 
dimension 2, coexistence occurs for an open set of values (Ai, A2) in R 2 . We conjecture that such 
a property holds in any dimension but our proof heavily relies on results of Durrett and Neuhauser 
(1991) that have been proved in d = 2 only. Before going into the details of the proof, we start 
by describing the setting of their results. First of all, let A 2 = 0. We recall that in this context, 
can be interpreted as a living tree, 1 as a burning tree and 2 as a burnt site. See Sect. 1 for more 
details. We start the process with one burning tree in a virgin forest and wonder whether the fire 
will go out or not. To do this, let B = (-L, L) 2 and, for any to € Z, B m = mLe\ + B, where 
ei = (1, 0) is the first unit vector. Then, denote by £ the subset of Z 2 given by 

£ = { (m, n) <E Z 2 , to + n is even } 

and say that a site (to, n) <G £ is occupied if the following two conditions are satisfied. 

(1) There are more than L 1/2 burning trees in B m at some time t€ [n T L, (n + 1) T L]. 

(2) There is at least one burning tree in B m at all times t £ [(n + 1) T L, (n + 2) T L]. 

Here F is a positive constant that will be fixed later. The following lemma is a comparison result 
of Durrett and Neuhauser (1991). It guarantees, in the setting of the forest fire model, that 
coexistence occurs for Ai greater than some critical value a c e (0, 00). 

Lemma 6.1 If A2 = and Ai > a c then T and L can be chosen so that the set of occupied sites 
dominates the set of wet sites in a one-dependent oriented percolation process on £ with parameter 

p=1 _g-36 

The next step is to show that the conclusion of Lemma 6.1 holds for small values of the parameter 
A 2 . To do this, take a site (to, n) £ £, and denote by B^ n the space-time region 

B^ n = B m x [nT L, (n + 2)T L] U B m+1 x [(n + 1) F L, (n + 3) L L}. 

Let H^ n be the event that the green particles do not give birth in the box B^ n . A straightforward 
calculation implies that the volume of is bounded by 14 T L 3 so that 

P(H^n) > e- 1AX - TL3 > I-6- 36 
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for A 2 > small enough. Denoting by B m n the set 

B m ,n = B m x [nTL,{n + 2)TL] U B m ^ x [{n + 1) T L, (n + 3) T L] 

and by i?^ „ the event that the green particles do not give birth in B mnl the same arguments 
obviously imply that 

P (#",„) > 1-6- 36 . 

On the other hand, Lemma 6.1 states that if (m, n) is occupied and H mn (resp. H mn ) occurs, 
the probability that (m + 1, n + 1) (resp. (m — 1, n + 1)) is occupied can be bounded from below 
by 1 — 6~ 36 . Putting things together, we can finally maintain that, for A 2 > small enough, the 
conclusion of 6.1 holds for p = 1 — 2 x 6 ~ 36 . 

To conclude the proof of Theorem 3, observe first of all that classical results about percolation 
imply that for such a value of p there exists a.s. an infinite cluster of wet sites in £. See for instance 
Durrett (1984), Sect. 10. Next, start the process from an initial configuration that has at least 
L 1/2 red particles in each box B m and extract a convergent subsequence of the Cesaro averages. 
The limit p we then obtain is known to be a stationary distribution. Moreover, by comparison 
with the percolation process defined above, fj, concentrates on configurations with infinitely many 
red particles. This concludes the proof of Theorem 3. 
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